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The Landau gauge quark propagator in momentum space is investigated using the C(a)-improved 
Sheikholeslami-Wohlert (SW) quark action with a tree-level mean-field improved coefficient c sw . We 
have studied the unimproved definition of the quark propagator, as well as two different tree- level 
O (a) -improved propagators. The ultraviolet behavior of the free lattice propagator is studied for 
each of these in order to establish which of them provides the most reliable description of the 
quark propagator up to the medium momentum regime. A general method of tree- level correction 
■ is introduced. This exploits asymptotic freedom and removes much of the trivial lattice artifacts 

' at medium to high momenta. We obtain results for the quark propagator which are qualitatively 

, similar to those typically used in quark models. A simple extrapolation of the infrared quark mass 

M(p 2 = 0) to the chiral limit gives 298 ± 8 ± 30 MeV, which is consistent with phenomenological 
expectations. 



oo 



C3 



PACS numbers: 12.38.Gc,11.15.Ha,12.38.Aw,14.70.Dj 



in 

> 
oo 

O 

r- . 

O I. INTRODUCTION 
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The quark propagator is one of the fundamental quantities in QCD. By studying the momentum-dependent quark 
mass, obtained from the scalar part of the inverse quark propagator, we can gain valuable insight into the mechanism 
of chiral symmetry breaking and its momentum dependence. The quark propagator is also used extensively as an 
input in Dyson-Schwinger Q| based model calculations of hadronic matrix elements |2|, || . Hence a lattice calculation 
of the quark propagator would enable us to check the validity of the models used in these calulations. There have been 
several recent studies of the quark propagator on a finer lattice with the aim of obtaining the light quark masses and 
. £^ | renormalization constants . Here we will focus more on the infrared and medium-momentum regime and extend 
some earlier preliminary work ||. For comparison with the present studies, some results for the quark mass function 
?-h ' using Kogut-Susskind fermions have recently been reported Q . 

The study of the quark propagator on the lattice is complicated by the explicit chiral symmetry breaking in the Wilson 
fermion action, and also by finite lattice spacing effects, which are large compared to those in the pure gauge sector 
0. For the gluon sector, on the contrary, one can achieve reliable results even with very coarse lattices using 0(a 2 )- 
improved actions together with mean- field improvement. These studies have shown, for example, that in Landau 
gauge the gluon propagator is enhanced at intermediate momenta and suppressed in the infrared to the point where 
it is almost certainly infrared finite [|| . 

Perturbation theory in a covariant gauge has a gauge-fixing parameter £, which corresponds to the width of the 
Gaussian average over the auxiliary field c(x) in the gauge fixing condition d^A^ix) = c{x). The choice of Landau 
gauge, i.e., £ = 0, corresponds to the zero width case, i.e., d^A^{x) = 0, which is the Lorentz gauge-fixing condition. 
Hence "covariant gauges" are actually Gaussian weighted averages over generalizations of the Lorentz gauge fixing 
condition. On the lattice, Landau gauge means that we have imposed the Lorentz gauge condition by finding a local 
minimum of the appropriate gauge-fixing functional [J7J. As for the previously cited studies of the gluon propagator, 
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the work reported here is done in the quenched approximation and without attempting to avoid Gribov copies, e.g., 
without attempting to project onto the fundamental modular region. Of course, in our finite ensemble of gauge field 
configurations no two Landau-gauge configurations will ever be Gribov copies of each other. However, the Landau 
gauge configurations will not be samples from a single connected manifold such as the fundamental modular region. 
This is an interesting area for future study. 

In a covariant gauge in the continuum the renormalized Euclidean space quark propagator must have the form 

a, \ = Z (^P 2 ) 1 m 

where we see that Z(n;p 2 ) = 1/A(n;p 2 ) and M(p 2 ) = B(n;p 2 )/A(fi;p 2 ). The renormalization point is denoted by /i 
and since we are interested in defining nonperturbative renormalization we use the standard momentum subtraction 
scheme (MOM) , which has the renormalization point boundary conditions 

Z(fi,fi 2 ) = 1 and M{pi 2 ) = m(/x) . (2) 

At sufficiently large /i in an asymptotically free theory like QCD the effects of dynamical chiral symetry breaking 
become small and to(/k) becomes the usual explicit chiral symmetry breaking running quark mass. At large Eu- 
clidean momentum scales (i.e., large fj,) the procedure for relating the parameters of the MOM scheme to the popular 
perturbative renormalization schemes (i.e., MS or MS) is well known. 

The renormalizability of QCD implies that the bare propagator is related to the renormalized one through the quark 
wavefunction renormalization constant Z 2 : 

S h ^(a;p) = Z 2 ( f i-a)S( f i;p) , (3) 

where a denotes the regularization parameter in some regularization scheme (such as the lattice or dimensional 
regularization). In a renormalizable theory the renormalized quantities become independent of the regularization 
parameter in the limit that it is removed (i.e., a — > on the lattice or e — » in a dimensional regularization 
scheme), while holding the renomalization point boundary conditions fixed in Eq. (^]). It immediately follows from 
the renormalization point independence of the LHS of Eq. (||) that for sufficently small a (i.e., in the scaling region) 
we have 

Z2(^a) Z(fi',p 2 ) *,f m 2\-M/ 2x A/U I 2n f.s 

Z 2 (/i';a) Z(p,p 2 ) 

for all p 2 . Hence the mass function must be renormalization point independent and a change of renormalization 
point is just an overall rescaling of Z(ji\p 2 ) by a momentum-independent constant, i.e., the LHS of the first equality 
in Eq. (|J). Hence, once the momentum-dependent renormalized propagator is known at one /i for all p, then it is 
immediately known for all [i. We can evaluate the constant needed to rescale Z(/j,;p 2 ) to Z(fj,';p 2 ) by evaluating 
Eq. (|) at p 2 = fi' 2 and using Z(ji'; ^' 2 ) = 1 [i.e., Eq. (§)] to give 



Z 2 (n;a) 



Z(^' 2 ). (5) 



Perturbative QCD chooses a renormalization scale /i close to the momentum scale characterizing the particular process 
of interest, e.g., /i 2 ~ Q 2 in deep inelastic scattering. This choice is made to ensure that perturbation theory will 
converge as rapidly as possible for the process of interest. The quark propagator used in such calculations is then 
S(fi;p) for p 2 near fj 2 , i.e., S peItulh {fJ,;p) = $ + m(/x)]. 

The tree- level quark propagator is the bare (i.e., regularized) quark propagator in the absence of interactions, i.e., 

S(°\a;p)= 1 , (6) 
i p + m u (a) 

where m°(a) is the bare quark mass. When the interactions with the gluon field are turned on then 

S^ip) - 5 barc (a;p) = Z 2 ( W o)S( W p) • (7) 

So we see that in the scaling region (i.e., for sufficiently small a) the measure of nonperturbative physics is the deviation 
of Z(/j,;p 2 ) from 1 and the difference of M(p 2 ) from the renormalized quark mass m(/z). As already mentioned, for 
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sufficiently large [i, becomes the running mass at the renormalization point fi, which is the basis of the studies 

performed in Rcf. Q. 

The purpose of the work reported here is to extract the full Z(fj,;p 2 ) and M(p 2 ) directly from a lattice calculation 
of the bare quark propagator S haTe (a; p) . It is of course 5 bare (a;p) that is calculated on the lattice. In reality we do 
not have the convenience of having an arbitrarily small a, rather we are faced with a lattice spacing which introduces 
lattice artifacts at medium to high momenta. In order to simplify the presentation of the data we will not explicitly 
introduce a renormalization point. Rather we will introduce for convenience the rcnormalization-point-independcnt 
combination 



(8) 



The regularization parameter dependence (i.e., the a-dependence) of Z(p 2 ) is not indicated for brevity, but is to be 
understood. We will develop a procedure for tree-level correction of the lattice artifacts in order to minimize their 
effect. 

The structure of the remainder of this paper is as follows: In Sec. |n] we describe the various 0(a) improved quark 
actions and propagators that we will study. In Sec. Ill we introduce our notation for the propa gato r in momentum 
space and derive the tree-level expressions appropriate for the actions that we consider. In Sec. IV we describe our 
methods to mi nimiz e lattice artifacts, including our tree-level correction scheme. Section ^ contains our numerical 
results: In Sec. V A we pres ent t he data f or th e pr opag ator without the tree-level correction; the effect of the tree-level 
correction is shown in Sec. VB; in Sees. VC and VD we present the results of fits t o a model function for the mass 



function M and extract the dynamically generated infrared quark mass; and in Sec. VE we discuss the possibility of 
finite volume effects on Z(p 2 ). Finally, in Sec. VI we present our conclusions and suggestions for further work. 



II. IMPROVED QUARK PROPAGATORS 



A systematic program of improvement || proceeds by adding all possible higher-dimensional local operators to 
the Lagrangian. When applied to the fermionic part of the QCD action, adding all possible gauge invariant local 
dimension-5 operators yields the following Lagrangian jl0[ [ll] : 

C(x) = C w - jCswaipa^F^if) + tr(F^F^) - b m am 2 ipip 

+c 1 mpp 2 i} + c 2 amAppij) . (9) 

Here for notational brevity we introduce the simple notation to for the lattice bare mass, i.e., m = m°(a). In this 
equation we have used C w for the standard Wilson Lagrangian density and (i/^c^aifxr^vF^i/; is the so-called 
"clover" improvement term. The sum of these two terms is often referred to as the Sheikholeslami-Wohlert (SW) 
action. That the action given by Eq. (^) is sufficient to remove all 0{a) errors has only been rigorously demonstrated 
for on-shell quantities. For gauge dependent quantities, it is an open question whether further, gauge non-invariant 
(but BRST invariant) terms must be added. We will assume that any such terms will be small. We shall follow the 
procedure used in studies of the gluon propagator ||, where it was seen that a combination of improved actions and 
tree-level correction gave reliable outcomes even at medium to high momenta. 

Since the Wilson action explicitly breaks chiral symmetry, the lattice bare mass should be taken as the so-called 
subtracted bare mass @ 

m= m°(a) = m - m c = - ( — ) , (10) 

a \2k 2k c J 

where too = (l/2na) — 4/a is the bare quark mass appearing in the Wilson action. At tree level, where interactions 
are absent, the quark condensate will vanish when the bare mass appearing in the action vanishes, i.e., when too = 
or equivalently when k = 1/8. In the interacting theory, n c is defined as the value of k at which the pion mass 
vanishes and m c = \/{2n c a) — 4/a is a nonperturbative fine-tuning correction needed to ensure that the bare mass 
to vanishes when the pion mass vanishes. The b g and b m terms correspond to a (mass-dependent) rescaling of the 
coupling constant and the mass respectively. Since we we will work in the quenched approximation we can set b g = 0. 
The parameter b m will be absorbed into a redefinition of the bare mass m here and we will comment on this later. 
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At tree level, the c\ and C2 terms can be eliminated by the following transformation of the fermion field jl2| , 

\j) ->V'= (l + b q am)(l - c q ap)ip 

i> -> i>' = (f + 6,om)^(l + c q a1jj) (ff) 

In general, beyond tree level, the improvement in the action must be combined with a corresponding improvement in 
the fermion field jll| , 

tP' = {l + b' q am)(l-c' q ap)TP + c n flip, (12) 

where the gauge dependent coefficient c„ is needed when we compute gauge dependent quantities, like the quark 
propagator. By choosing the correct improvement coefficients for the field, the c\ and C2 terms may again be eliminated. 
We note in passing that the coefficients b' q and c' q were recently calculated at one-loop level |13[, while c„ is still 
unknown. However, here we will be restricting ourselves to tree-level 0(a)-improvement throughout for the coefficients 
b' q , c' q , and c' n . In that case, the 0(a) improved action and fields after the transformation in Eq. ( p"l| ) have b' q — > b q = j 
and c' q — > c q = \, and c„ = 0. We will use the tree- level mean-field improved value for c sw and the nonperturbatively 
determined value of k c to determine the lattice bare mass m in terms of k. Although we will use tree-level improvement 
formulae for our quark actions and propagators, it is more appropriate to use m than m,Q in these. Although there is 
an apparent inconsistency in using tree- level values for b' q and c' and the mean-field improved value for c sw , we have 
numerically verified that using mean-field improved values for b q and c' q makes no significant difference in practice. 

The tree- level C(a)-improved propagator can then be defined as 

S(x, y) = (V'(x)^'(y)) = ((1 + b q amf{\ - c q ap(x))S (x, y; U)(l + c q af(y))) , (13) 

where b q — c q ~ 1/4 and where Sq(x, y; U) for a given configuration U is simply defined as the inverse of the fermion 
matrix, 

M (x; U) =p w {x; U) + {ia/^c^a^F^ix) + m =p{x; U) + m + 0{a) , (14) 

where p\y{%] U) is the lattice Wilson-Dirac operator and-p(x; U) is the usual continuum covariant derivative. There- 
fore, So(x,y; U) will always satisfy the relations 

[p(x:U)+m]S (x,y;U) = 8{x - y) + 0{a) 

S (x,y;U)[- < p(y;U) + m] = 6(x ~ y) + O(a) (15) 

The "unimproved" quark propagator Sq will be defined here to be that arising from the SW action consisting of the 
Wilson term and the clover term, but with no other corrections. Hence, Sq is then given by the ensemble average of 
S Q (x,y; U) 

S (x,y) = (S (x,y;U)) . (16) 



We will denote the tree-level C(a)-improved quark propagator obtained from Eq. (13) as the improved "rotated" 
propagator Sn(x,y), which is 



am, 



2 



a 



S R (x,y) = {S R (x,y;U)) = {(l + —) l--Mx) S (x,y;U) l + -p(y)). (17) 
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(18) 



We can use Eq. (|15|) to obtain another, simpler expression for the improved propagator from Eq. (fL3| 

S(x, y) = (l + 2(b q + c q )am)S Q (x,y) - 2ac q S(x - y) + 0(a 2 ) 
= (1 + am)S (x -y)- ^6(x - y) + 0(a 2 ) , 

where we have used the fact that b q = c q = 1/4 here. We define the corresponding version of the tree-level O(a)- 
improved propagator as 

Si(x -y) = (l + am)S (x - y) - -8{x - y) (19) 

If we are only interested in on-shell improvement, e.g., hadronic matrix elements, the 5-function can be ignored. 
However, it is essential if we are considering off-shell properties such as the quark propagator in momentum space. 

In summary, we see that both Sr and Si are tree- level improved definitions of the SW-clover (i.e., the 'Wilson plus 
clover') propagator Sq. However Sr and Sj will have different 0(a 2 ) errors in general. This will become an important 
consideration when we later attempt to minimize lattice artifacts. 
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III. MOMENTUM SPACE PROPAGATOR 



The momentum space quark propagator is given by 

S(p) = J2< ripX S( x ,°) (2°) 

X 

As is appropriate for fermions we will be using periodic boundary conditions in the spatial directions and antipcriodic 
boundary conditions in the time direction. Hence, the available momentum values for an Nf x N t lattice (with iVj, N t 
even numbers and i = x,y, z) are 

Pi = T7 — ( n i 7T ) 5 m =1,2,- •• ,Ni (21) 



N t a V 2 

We will also for notational convenience define the following 'lattice momenta': 

fc p = - sin(p M a) (23) 

2 / 

L = - sin(p„a/2) = J 1 - cos(p„a) (24) 

a a V 

which differ by 

a 2 A^P-fc 2 ^5>4 +0(a4) (25) 

In the continuum, the quark propagator has the general from given by Eq. (Q). On the lattice it is convenient to work 
with the dimensionless quark propagator S(p) = 5 bare (a; p) / a. We expect the lattice bare quark propagator to have 
a similar form to its continuum equivalent, but with the lattice momentum ft replacing p 1 , which can be appreciated 
by referring to the tree level lattice propagators to be given later. Because of hypercubic lattice artifacts Z L and 
the dimensionless M L will be functions of p M rather than p 2 . Hence we have for the dimensionless lattice bare quark 
propagator the form 

o (] _ 1 _ Z L (p) _ L —ia ft + Al L (p) 

[P! ia ftA{p) + B{p) ~ iaft + M L (p) ~ {P> a 2 k 2 + (M L ) 2 (p) [ ' 

In the limit a — ► the continuum form will be recovered. The (dimensionless) lattice functions A(p) and B{p) can 
then easily be extracted from the inverse dimensionless lattice quark propagator, 

aM-f^-JS^'M) (28 » 

In practice, however, it is easier to extract these functions without inverting the propagator. It is easily verified that 

A{p) = k 2 a 2 A 2 {p)\ B 2 (p) (29) 
B ^ = k 2 a 2 A 2 (p) + B 2 (p) (30) 



where we have defined 



^mW^WW B(p) S ±-tr S(p) (31) 
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A. Tree-level expressions 



As was done in the introduction in Sec. I we will use the superscript (0) to denote the tree-level versions of each of 
the quark propagator definitions and actions considered. These are simply the propagators that would be obtained 
from the various definitions when the interaction with the gluons is turned off, i.e., when all the gluon links are taken 
to be the identity. We will also be writing m rather than mo throughout, since at tree level the two are identical. The 
dimensionless SW fermion propagator at tree level is identical to the pure Wilson propagator and is given by (lj, [ll| 



Wa 2 



k 2 a 2 



fk 2 a 2 



(32) 



The tree-level form of the tree-level 0(a)-improved propagator Sj is given by 



(33) 



If we write 



we find 



(sf (?)) 1 =i}6aA < ?\p) + B?\p), 



4%) - -jjrtr 



(sPiP))' 



/k 2 a 2 = l + 0{a 2 ) 



Bf\ P ) 



. , ma 2 

///-/ 1 ; 1 - — + 0{a 2 



(34) 

(35) 
(36) 



We see that the quark mass gets an 0{a) correction. The purpose of the improvement term b m in the action (Q) was 
to cancel this change in the bare mass m. However, by omitting this correction we have have simply absorbed it into 
a redefinition of m. A^°\p) is equal to unity up to C(a 2 ), as expected. The details of the derivation are given in the 
appendix. 

It is also useful to write the propagator in the following way, 



1 



Z(°)(p) 



i }ta + ma + aAM (0) (p) 



(37) 



The analytic expressions for Z^(p) and AM^(p) for both the improved propagators Si and Sr are given in the 
appendix. To illustrate the behavior of these tree-level functions we show in Fig. |l| the forms of and AM' ' 

for both of our improved actions Si and Sr. The horizontal axis is pa = ^Jp 2 + p 2 + p 2 , + p 2 , where the possible 

values of the momenta are those given in Eq. (p2]). Since Z^(p) deviates from 1 and AM^°\p) deviates from zero 
at medium to high momenta, it is immediately obvious that the finite a effects are very large and we will need some 
method for taking care of them if we are to obtain physically meaningful results. The tree-level behavior is particularly 
pathological for Si(p), with finitc-a effects of several hundred percent appearing in Z(°), and AM' ' being many times 
larger than m, and negative. The spread in the points is due to hypercubic artifacts, since on the lattice Z (0) and 
AAf(°) are functions of p^ and not p 2 . The finite-a effects in Sr are much more mild and offer the hope that they 
might be partially compensated for. Clearly in the limit a — > we recover the continuum result where Z^°'{p) = 1 
and AM(°)(p) = for all p. 



IV. ANALYSIS 



A. Tree- level correction 



Recall that the quark propagators calculated on the lattice are actually the bare quark propagators, which become the 
tree-level propagators when the interactions are switched off. We know that QCD is asymptotically free, which means 
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FIG. 1: Plots of the analytic functions (top) and oAM* 9 ' (bottom) versus the mom e ntum pa for Si{p) (left) and Sr(p) 
(right). The results shown here are obtained from the analytical expressions in Eqs. ( A8)-(<\15). 



that at sufficiently high momentum values the bare quark propagator should approach the tree-level quark propagator 
up to logarithmic corrections, i.e., on the lattice for large momenta we should find S(p) — > S^(p) up to logarithms. 
The deviation of these from each other is a direct measure of the nonperturbative effects due to the interactions felt 
by the quarks. Hence, we are here primarily interested in studying the deviation of the quark propagator from its 
tree-level form. 



We will attempt to separate out the tree-level behavior by writing 

1 



S-\pa) 



ia f6 + aM(pa) + aAAf (0) (pa) 



(38) 



If asymptotic freedom holds for the momentum range we are considering, we should expect that Z(pa) — > 1 and 
M(pa) — > m (up to logarithmic corrections) for large p. 

Equation ( |38| ) can be rewritten to yield expressions for Z(pa) and M(pa) in terms of the functions Z L and M L (or, 
equivalently, A and B) defined in Eq. (p6[), 

(39) 
(40) 



Z L (pa) 

Z{pa) = zMh 

aM(pa) = M L (pa)~ aAM {0) (pa) 
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We refer to the functions Z and M obtained in this way as the tree-level corrected forms of the lattice quantities Z L 
and M L . 

It is important that one not be confused by the different uses of the expression "tree-level". Firstly, there is an 0(a)- 
improvement which was only implemented at tree-level rather than having the improvement coefficients determined in 
some nonperturbative way. Secondly there are tree-level propagators which are the bare propagators when there are no 
interactions. Thirdly, we have just now introduced our method of tree-level correction, which will hopefully minimize 
the finite-a errors in our extraction of the quark propagator from the lattice. Because the tree-level behavior of Si, 
i.e., Sj°\ is much worse at medium and high momenta than that of Sr, we anticipate that our tree-level correction 
may not be adequate in that case. We therefore expect the tree-level correction method to work significantly better 
for Sr than for Sj. 



B. Cuts 



Even after the tree-level correction has been performed, there will still be anisotropies in the data, resulting from 
finite a effects beyond tree level. To remove these, we select momenta lying close to the diagonal in momentum space. 
We choose the diagonal because finite-a hypercubic artifacts will be minimized at a given pa when the momentum 
is approximately equally spread among the four momentum components. Ideally, one should attempt an a — ► 
extrapolation, but given the available data we will see that this cut on the data removes most of these artifacts. We 
define the distance of a point from the diagonal by 

Ap=\p\sm9(p) 1 (41) 

where the angle 9(p) is given by 

cos%) = ^, (42) 
\P\ 

and h = i(l, 1, 1, 1) is the unit vector along the diagonal. We select momentum values such that Ap < 7r/8a, ie. 
within one unit of spatial momentum from the diagonal. We refer to this selection as the "cylinder cut", since the 
momenta selected lie within a cylinder around the diagonal in momentum space. 



V. RESULTS 



The quark propagator is calculated at j3 = 6.0 on a 16 3 x 48 lattice, using the tree-level mean-field ('tadpole') improved 
value c sw = 1.479. For this action with these parameter values, k c is found to be k c — 0.1392 [fill. Two values for 
k were used: k — 0.137, corresponding to ma — 0.0603, and k — 0.1381, corresponding to ma — 0.031. The lattice 
spacing as determined from string tension measurements in the gluon sector at (3 = 6.0 is a = 0.106 ± 0.002 fm or 
equivalently 1/a = 1.855 ± 38 GeV and so the values for the quark masses are m = 112 MeV and m = 57.5 MeV 
respectively. 

The configurations were fixed to Landau gauge with an accuracy of 8 = Y^, x ,u I^VA^ 2 -)! 2 < 10 -12 . At k — 0.137, 
we have generated both So and Sr. At n = 0.1381, only So was generated. Si is easily constructed from So- We 
have used the tree-level values for the coefficients b q and c q as previously stated, rather than the mean-field improved 
values, as the difference between the two is negligible compared to the 0(a 2 ) and higher effects which the tree-level 
correction scheme attempts to minimize. We have explicitly verified that replacing the tree-level values for b q and c q 
with the mean-field improved values makes only negligible difference. However, in any future study it would clearly 
be preferable to make consistent use of mean-field improved or non-perturbatively determined (as far as they are 
available) improvement coefficients throughout. All the results shown for Sr are for 20 configurations, while the 
results for Si are for 499 configurations, unless otherwise specified. 

As a further check on our results, we have also analyzed 60 configurations at j3 — 5.7 on a 12 3 x 24 lattice, for 
k = 0.13843 and k = 0.14077, corresponding to ma = 0.128 and 0.068 respectively. Here, k c = 0.1432. In this case, 
all three propagators were generated for all configurations. We will not explicitly show these results here but will 
comment on their relevance in our later discussion. 
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FIG. 2: Z L (p) — 1/A(p) as a function of momentum p for Si(p) (left) and for Sr(p) (right) with the bare quark mass 
corresponding to k = 0.137. No tree-level correction has been made and no data cuts have been applied. 




FIG. 3: M L (p) for Si(p) (left) and for Sr(p) (right) with the bare quark mass corresponding to c = 0.137. No tree-level 
correction has been made and no data cuts have been applied. 



A. Uncorrected data 



Let us first see what happens when we use the naive formulae for Z and M without implementing the tree-level 
correction, i.e., we first consider Z L and M L . Since Eq. ( [l5| ) holds precisely configuration by configuration, Eq. ( |Tg| ) 
should be satisfied non-perturbatively. However, the 0(a 2 )-term can be quite large. In Fig. ||we show Z L (p) — 1/A(p) 
as a function of pa using Si(p) and Sr(j>) respectively, while in Fig. |^ we present M L (p) for the same two cases. 
Comparing these figures with the tree-level behavior shown in Fig. [j] we see that finite-a errors completely dominate 
these (uncorrected) quark propagators at medium and high momenta. Only in the infrared, below pa < 0.8, might 
we be able to extract physically significant information. 

A comparison with the 'unimproved' SW propagator So shows that both Si(p) and Sr(p) are considerably better 
behaved in the infrared than the 'naive' propagator So(p). In particular, the mass function is a decreasing function of 
pa up to pa ~ 1, which is what one would expect from asymptotic freedom. This is not the case for Sq(p), which begins 
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FIG. 4: Z(p) (left) and aM(p) (right) for our preferred form of the improved propagator Sr(p) at k — 0.137. The lattice data 
shown are obtained using the tree-level correction defined in Eq. (^i|) but without any cuts. 



to increase monotonically in the infrared at pa ~ 0.4 as does its tree-level form. Also, the values for Z L (p) = 1/A(p) 
and M L (p) agree for 5^ and Si within errors up to pa ~ 0.8, while the values obtained from So(p) are significantly 
different from the improved values even at low momenta. 

At momenta above pa ~ 1 the 0(p 2 a 2 ) and higher terms dominate and it is impossible to extract any meaningful 
information from these uncorrected data. This can be appreciated most dramatically by the way the lattice data 
diverge with increasing momentum for the two improved (but uncorrected) propagators. 



B. Tree-level corrected data 



Since, as we saw in the previous section, the tree-level form completely dominates the high-momentum data, we may 
hope that by factoring out this behavior we will get something which lies close to the continuum asymptotic form. In 
the low and intermediate momentum region we may then be able to extract the physical, nonperturbative behavior 
of the functions Z(p) and M(p). 

When applying this correction, we find that there is a dramatic improvement in the behavior with pa of all three of 
our forms of the propagator, i.e., for Sq, Si, and Sr. However, the pathological behavior of Sj°\p) at high momenta 
gives rise to a cancellation of large terms in the subtracted mass, leading to a behavior for the mass function which 
is clearly at odds with the expectation from asymptotic freedom. Thus, as expected, the finite-a errors in Si are 
simply too large to be corrected by our simple tree-level correction procedure. As previously noted, our unimproved 
propagator Sq behaves poorly even at very low momenta and cannot therefore be trusted. It is therefore desirable to 



use the definition Sr for the improved propagator) 19 and to apply our tree-level correction to that. The results for 
Z(p) and M(p) for our preferred propagator Sr(p) are shown in Fig. || as functions of pa for k = 0.137. We see that 
the medium to large momentum behavior has been dramatically improved by our tree-level correction procedure as 
expected, i.e., it behaves in a way reasonably consistent with the expectations of asymptotic freedom. The spread of 
the lattice data due to hypercubic artifacts is somewhat reduced but has not been eliminated. 

In Fig. ^ we show the lattice results for Z and M for all three definitions of our quark propagator, after implementing 



both the tree-level correction and the cylinder cuts described in section IV B| . The tree-level correction is clearly 



failing for Si as is evident from the behavior of the mass function M(p). Although the behavior of the unimproved 
propagator So has been considerably improved, as we previously observed it cannot be trusted even at relatively 
low momenta and so it must be discarded. The apparent difference in the behavior of Z between the two improved 
definitions of the propagator Si and Sr, even at relatively low momenta, is at first sight puzzling. However, it 
must be recalled that Z(p) = Z2((J.;a)Z([i;p) and that it is actually Z(p,;p) that we should be comparing for the 
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FIG. 5: The tree-le vel co rrected functions Z(p) (left) and aM(p) (right) for all three propagators, after performing the cuts 
described in section IV B . The values for Sr are obtained from 20 configurations; for the two other propagators the data 
are from 300 configurations. Note that the Z(p) functions need to be scaled to agree at some "safe" momentum scale (e.g., 
fia ~ 0.4) before being compared. 



different actions. Different actions will in general have different values of the renormalization constant Z<x{p,\ a). If we 
renormalize at some "safe" momentum scale where we would expect both improved propagators to be reliable, e.g., 
fia ~ 0.4, then the apparent difference is much reduced except at medium to high momenta where Si can no longer 
be trusted. Analysis of the data at j3 = 5.7, which corresponds to a coarser lattice, and hence larger finite-a effects, is 
also consistent with this interpretation. Below pa ~ 0.6, the values for M{p) agree within errors for the two versions 
of the improved propagator. In particular, the value for the infrared mass M(p — > 0) agrees well. In contrast, and 
not surprisingly, the unimproved propagator Sq yields a mass which is 3-4<r higher. We again clearly see from this 
plot the poor high- momentum behavior of M(p) arising from the inexact cancellation of large finite-a errors in Si. 

Fig. |^ shows M(p) calculated from Si for the two values of the quark mass. In the infrared region, the mass changes 
only slightly as the bare quark mass is halved, pointing to a dynamically generated 'constituent' quark mass in the 
chiral limit, which we will later estimate. The function Z(p) was found to be insensitive to the quark mass. 



C. Model fits 



In order to try to parametrize its behavior, the mass function M{p) was fitted to the simple model analytical form 

aM(pa) = - A2 + m uv , (43) 

where k was defined in Eq. We have fitted to data in the window < ka < P, with P varying between 0.7 and 

1.4, in order to verify that the parameters are insensitive to the fitting window P. Since Sr is far better behaved 
than Si at higher momenta, all the fits have been performed to the mass function M(p) extracted from our preferred 
propagator Sr. 

The parameter values for k = 0.137 are shown in table |. All the fits give a value for m uv which is consistent with 
0. This is due to the fact that we have not completely removed these lattice artifacts from M(p) at intermediate and 
large momenta. This indicates that at /? = 6.0 with our preferred improved action and propagator we still do not 
have sufficient control of ultraviolet lattice artefacts that would allow us to extract the ultraviolet running mass [H . 
Combining the fit parameters for all the fits gives a value for the infrared quark mass of aM- u = 0.211 ± 0.008. 
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FIG. 6: The tree-level corre cted m ass function M(p) from Si, for k = 0.137 (circles) and k — 0.381 (crosses), after performing 
the cuts described in section IV B. The increase in M(p) for pa > 1 is an indication of the difficulty of accurately subtracting 
off the tree-level mass function for this definition of the improved propagator. 57 has been used in this case because we lack 
data for more than one quark mass for our preferred definition Sr; however, the mass functions agree well in the infrared. The 
data shown are from 300 configurations. 



TABLE I: Parameter values for best fits to the form of Eq. (|l3j) in the window < ka < ap max , for different values of p max , at 
K = 0.137. 
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D. Infrared quark mass 



A quantity of intrinsic interest is the mass function M(p) at zero momentum in the 'chiral' limit m — > 0. This gives a 
measure of the dynamical chiral symmetry breaking in the system, and is related to the order parameter of dynamical 
chiral symmetry breaking, the chiral condensate (ipip), as well as to the concept of the 'constituent quark mass' used 
as input in various quark models. 

Since we have only computed Sr for one value of the quark mass, we must use the data from Si to perform the 
extrapolation to m = 0. Recall that at low momenta the two actions give consistent results for M(p ) within errors. 
The results are shown in Fig. |L The data point from Sr at k = 0.137 is also shown, giving an indication of systematic 
uncertainties. We find a value for M(0; m = 0) of 298±8±30 MeV, where the second set of errors is an estimate of the 
systematic uncertainty coming from the difference between Sr and Si. This is reasonably consistent with commonly 
used values for the constituent quark mass. 



E. Finite volume effects 



To determine whether the infrared suppression of Z{jp) is real or a finite volume effect, we can look for anisotropy 
in the infrared. Since the temporal extension of the lattice is three times the spatial extension, the finite volume 
will affect spatial momenta differently from timelike momenta, giving an indication of the size of (anisotropic) finite 
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FIG. 7: The lattice results for Z(p) = Z^iji; a)Z(n;p) and M(p) for our preferred form of the quark propagator Sr(p), after 
both the tree-level correction and the cylinder cut. The vertical scale for Z([i;p) is determined from the above by dividing it by 
the necessary renormalization constant (i.e., ^2(^1 °)) to ensure that Z(fi; (j 2 ) — 1. These are the central results of the studies 
reported here. The bare quark mass used here was m = 112 MeV and hence we conclude that M(p) is not reliable at momenta 
above approximately 1.5 GeV. 
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FIG. 8: The infrared value of the quark mass function Mi r = M(p = 0), obtained by extrapolating M(p) to pa = for two 
different bare quark masses, m = 57.5 and 112 MeV. The crosses denote the two values obtained from Si from 300 configurations 
and two bare quark masses, while the square is the value obtained from Sr for a single quark mass. The burst indicates the 
chirally extrapolated value of Mi r obtained by a simple straight line interpolation. 



volume effects. 

Figure || shows the infrared behavior of Z(p), with momenta in different directions plotted separately. We see that the 
finite volume anisotropy, although not negligible, is not sufficient to explain the infrared suppression. This indicates 
that the suppression is either due to isotropic finite volume effects, or is a real physical phenomenon. In model Dyson- 
Schwinger equation studies (!]] the dynamically generated quark mass is typically associated with a dip in Z(p) similar 
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FIG. 9: Z(p) from from the propagator Si(p), from 70 configurations at k, = 0.137, i.e., m = 112 MeV. The squares denote 
purely timelike momenta, while the diamonds denote points with one unit of spatial momentum. The fancy squares are points 
with half a unit of timelike momentum (ie, nearly purely spatial momentum). These data appear to indicate that the infrared 
suppression of Z(p) is not a finite volume effect. 



to what is seen here. There is no discernible anisotropy in the data for the mass function M(p) at low momenta. We 
therefore conclude that finite volume effects for M(p) are almost certainly negligible. 



VI. CONCLUSION AND FURTHER WORK 



We have presented initial lattice results for the momentum-dependence of the quark propagator after implementing 
a tree-level correction procedure. At high momenta, quarks are asymptotically free and so the quark propagator 
approaches its tree-level behavior. We make use of this fact to subtract off and factor out the tree-level behavior, 
replacing it with what should be a more continuum-like medium and high momentum behavior of the quark propagator. 
This approach can only work reliably when the tree- level finite-a effects arc not too large, i.e., when the tree-level 
propagator corresponding to the action of interest is reasonably behaved at medium and high momenta. The tree-level 
correction was seen to dramatically improve the data for the preferred definition of the improved quark propagator Sr. 
The relatively poor behavior of the tree-level corrected Si is due to the large tree-level finite-a effects which require 
fine tuning to subtract off correctly. The unimproved propagator was seen to be unreliable even at low momenta 
(pa ~ 0.4) and so cannot be trusted even after tree-level correction. 

Although the ultraviolet behavior of the quark propagator is clearly improved, it remains an open question whether 
there exists a momentum window where the lattice data are reliable and perturbation theory is valid. One way of 
checking this would be to calculate the propagator at three or more different quark masses and attempt a chiral 
extrapolation of the mass function in the intermediate momentum region. If perturbation theory is valid in this 
region, the mass function should extrapolate to zero. 



Residual lattice artifacts may also be investigated by studying the chiral Ward identity JL7, 18 , which should be valid 
at all momenta. This involves computing the pseudoscalar vertex, and while it falls outside the scope of this initial 
study, it should be included in future studies of the quark propagator. In Ref. Jul the Ward identity has been studied 
for a slightly different action to ours, and verified for momenta up to pa < 1. 

The central results of this work are summarized in Figs. [7| |s| and ||. Fig. represents the best estimate from our 
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currently available data of the nonperturbative behavior of the quark propagator and is based on our preferred 
quark action corresponding to the Sr propagator. Fig. || is our extraction of the dynamically generated infrared or 
'constituent-like' quark mass. Finally Fig. | gives an indication of the magnitude of finite volume effects in Z(p) 
compared to the nonperturbative effects. 

For pa < 1, we see that M(p) falls off with p as expected. The values obtained from Sr and from Si are consistent, 
while those for the unimproved propagator So differ significantly. The infrared mass -M(O), which can be thought of 
as analogous to a 'constituent quark mass', appears to approach a value of 298 ± 8 ± 30 MeV in the chiral limit. 

We also find a significant dip in the value for Z(p) at low momenta. It must be remembered that the curves for 
Z(p 2 ) for different actions need to be rescaled to agree at some "safe" low momentum renormalization point before 
comparing them. The finite volume anisotropy is much smaller than the apparent infrared suppression. We cannot 
explicitly rule out large isotropic finite volume errors, although based on experience with earlier gluon propagator 
studies this seems unlikely. However, a larger volume is needed to completely resolve this issue. 

Since in this inital study, we have used the mean- field improved value for the clover coefficient c sw , and tree- level 
improvement for the fermion fields, the quark propagator still has some residual 0{a) errors as well as 0(a 2 ) and 
higher order errors. To remove the residual C(a)-errors it would be necessary to compute the non-perturbative values 
for the coefficients b' q ,c' q and c n . Repeating these calculations at a different lattice spacing and with other improved 
quark actions is also essential to get reliable results for the quark propagator and in particular for the quark mass 
function at medium to high momenta. These studies are currently underway. 
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APPENDIX A: TREE-LEVEL EXPRESSIONS — DETAILS 



The dimensionless Wilson fermion propagator at tree level is 



s?Hp) = 



-i jia + ma + ^k 2 a 2 



k 2 a 2 



ma + 7}k 2 a 2 



D 



= — —i tea + ma 



-k 2 a' 



(Al) 



Since the SW term is proportional to the gauge field tensor, it vanishes at tree level, so this expression also holds true 
for the SW action. The tree-level 'improved' propagator Si is given by 



Sf\p) = (l + ma)S^>(p) 



(0), 
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and the inverse propagator is 
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If we write 



we find 
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If we write the propagator according to Eq. (p 
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we find 
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where 



Bj = ma 



i 2 a 2 a 4 Ak 2 a 4 'k 4 



We have defined the rotated propagator Sjt(x,y) as 



S a (x,y) = (1 + ™) (l - ^(i)J S (z,y) (l + ^(y) 



At tree level, the fourier transform of this is 
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We can then write 
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j,From this we find the expressions for Z^ R = ^/A R and aAM^ = Z R B R — ma, via 
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Comparing these expressions with those of Eqs. 
are much smaller than for Sj. 



t) and (|Aq), we clearly see that the tree- level 0(a 2 ) errors in Sj 
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